The real shape of non-Gaussianities 

Antony Lewis^ 

^Department of Physics & Astronomy, University of Sussex, Brighton BNl 9QH, UIV] 

(Dated: March 15, 2012) 

I review what bispectra and trispectra look like in real space, in terms of the sign of particular 
shaped triangles and tetrahedrons. Having an equilateral density bispectrum of positive sign cor- 
responds to having concentrated overdensities surrounded by larger weaker underdensities. In 3D 
these are concentrated density filaments, as expected in large-scale structure. As the shape changes 
from equilateral to flattened the concentrated overdensities flatten into lines (3D planes). I then 
focus on squeezed bispectra, which can be thought of as correlations of changes in small-scale power 
with large-scale fields, and discuss the general non-perturbative form of the squeezed bispectrum 
^SJ and its angular dependence. A general trispectrum has tetrahedral form and I show examples of 

T-H what this can look like in real space. Squeezed trispectra are of particular interest and come in two 

f^ forms, corresponding to large-scale variance of small-scale power, and correlated modulations of a 

04 bispectrum. Flattened trispectra can be produced by line- like features in 2D, for example from cos- 

« I mic strings, and randomly located features also give a non-Gaussian signal. There are relationships 

Cu between the squeezed types of non-Gaussianity, and also a useful interpretation in terms of statis- 

^^ tical anisotropy. I discuss the various possible physical origins of cosmological non-Gaussianities, 

^^ both in terms of primordial perturbations and late-time dynamical and geometric effects. 



O I. INTRODUCTION 

u 

^ Most of the fields that we observe in cosmology are significantly non-Gaussian: the large-scale matter density is 

O, strongly non-Gaussian due to non-linear growth, and even the CMB is significantly non-Gaussian on small scales 

I due to non-linear effects. Current evidence suggests that primordial fluctuations at the beginning of the hot big 

2 bang were very close to Gaussian, but primordial non-Gaussianity is still a possibility and, if observed, would be a 

■^-^ way to rule out wide classes of early-universe models. Understanding non-Gaussianity is therefore crucial to extract 

jrt the most information from cosmological observations. Since the initial fluctuations were close to Gaussian, the non- 

I— ' Gaussianities can often be treated perturbatively, in which case the bispectrum and trispectrum contain most of the 

additional information (though at very late times where strongly non-linear processes are involved a deeper analysis 

may be required). Here I will focus on non-Gaussianity as described by a bispectrum and trispectrum. I will show 

what fields with various qualitatively different types of non-Gaussianity look like in real space, and give some general 

^^ results for the form of the squeezed bispectrum, its angular decomposition, relation to the trispectrum, and the 

^^ relation between anisotropy estimators and local non-Gaussianity. 

lO This paper does not aim to review the details of observational analysis or physical modelling of non-Gaussianities; 

r — there are many excellent reviews and references available for further technical and physical model details (for a small 

f^ sample see e.g. Refs. [TMlOj). I will mostly focus on generalities rather than specific models, though I mention specific 

T-H cases of cosmological interest as we go along. For simplicity I restrict to scalar fields in flat space; the generalization 

^~~^ to 2D fields on the sphere (e.g. the full-sky CMB), a non-flat background, or tensor flelds (e.g. galaxy lensing and 

^ CMB polarization) is conceptually simple, though can be technically significantly more complicated. 

• ^ I will assume a standard statistically isotropic and homogeneous background cosmology. In this case the wavevectors 

rN describing a particular non-Gaussian configuration of modes must sum to zero: for a bispectrum there are three modes 

^ with wavevectors that close to form a triangle; for a trispectrum there are four wavevectors that define four edges 

of a tetrahedron. By "shape" , in this paper I will mean a particular configuration of fixed-length wavevectors, e.g. 

a specific bispectrum triangle. A general bispectrum or trispectrum gives the full wavenumber and configuration 

dependence, i.e. the amount of signal expected in each possible shape. For example an 'orthogonal' bispectrum does 

not correspond to a particular triangle shape, but instead a particular wavenumber dependence of the signal in each 

shape. The more general wavevector dependence of the bispectrum and trispectrum has been discussed extensively 

elsewhere O O [H]; instead of showing the magnitude of the signal as a function of the wavevectors, I instead focus 

on what the signal looks like in real space for specific configurations. This can be a useful aid to understanding which 

kinds of physical process generate the different shapes. 
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II. GAUSSIANITY AND THE POWER SPECTRUM 

Before discussing non-Gaussianity, it is worth quickly remembering the key features of Gaussian fields. In particular 
we are usually interested is statistically isotropic and homogenous universe models, and hence in fields that have these 
symmetries. For simplicity I shall focus mainly on scalar fields 2D flat space, for example a slice through the matter 
density field or a small patch of the CMB, but almost everything generalizes to other cases such as full-sky observations. 
Assuming we can measure a field r(x) as a function of position, in fiat space these can be Fourier transformed and 
written as 

where N is the number of dimensions. Statistical homogeneity and isotropy means that the statistical properties of 
the field must be unchanged under translations and rotations r(x) -^ T{x'), so (r(x)T(x')) can only be a function 
of the invariant separation between the points |x — x'|. This implies that the covariance of the field is determined by 
a power spectrum depending only on fc = |k|: 

(T(ki)T(k2)) = <5(ki + k2)P{h). (2) 

For a small patch of the CMB, the power spectrum is just C/, where I — k. The delta function says that modes with 
different wavevectors are completely uncorrelated: knowing the sign of T(ki) tells you nothing about the likely sign 
of T(k2). From the power spectrum the only thing we know is the variance of each individual mode, which from the 
assumption of isotropy is the same independent of the orientation of the mode. 

A purely Gaussian statistically homogeneous and isotropic field is fully described statistically by its power spectrum. 
However more interesting fields are possible that are also statistically homogeneous and isotropic, with the non- 
Gaussian statistics described by a series of higher-point correlation functions. 

III. BISPECTRUM 

The first non-Gaussian signal to consider is a bispectrum, corresponding to a three-point correlation, or in Fourier 
space a correlation between three different mode wavevectors. We are still interested in statistically homogeneous and 
isotropic fields, which imphes the statistics are described by a reduced bispectrum b{ki,k2, k^) that depends only the 
lengths of the wavevectors: 

(T(ki)T(k2)T(k3)) = 7^4]^'5(ki + k2 + k3)6(fci, fc2, ks). (3) 

From now on when referring to the bispectrum I will mean the reduced bispectrum as defined here; an analogous 
definition applies on the full sky. The delta-function here means that the 3-mode correlation is zero unless the 
wavevectors sum to zero: they form a triangle. If there is a non-zero bispectrum, modes with different wavevectors 
are not independent: if we measure T(ki) and T{\s.2), the sign of the bispectrum b{ki,k2, k^) then tells us which sign 
of T(k3) is more likely. Positive sign gives positive skewness (tail of very high values), negative sign gives negative 
skewness (tail of very low values) . What this looks like in real space depends on the shape of the triangle (the relative 
lengths of the different wavevectors). 

A. Equilateral and flattened (folded) triangles 

The first case I consider is equilateral triangles, where the lengths of the three sides of the triangle are the same, 
ki = k2 = k^. If there's a non-zero equilateral bispectrum, what does the field look like in real space? To answer 
this we can consider taking the r(ki) and T{\^2) components of the field, and then ask what the bispectrum tells us 
about T(k3). Depending on the relative sign of T(k3), a field consisting of these three modes looks rather different 
— see Fig. fl] The sign of the bispectrum tells us which sign of r(k3) is more likely, in other words whether we are 
more likely to have small regions of concentrated overdensity [h > 0) or regions of concentrated underdensity (6 < 0). 
As can easily be imagined, such patterns can be obtained by locally moving matter around, for example concentrated 
overdensities can form by gravitational collapse, and thus equilateral non-Gaussianity is likely to be present in any 
field undergoing local non-linear dynamical processes. 

A bispectrum is determined by three wavevectors which always lie in a plane. In 3D, the modes we are considering 
correspond to plane waves, and the concentrated overdensities correspond to filaments. These are precisely what form 
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FIG. 1: Equilateral bispectrum: a field can be decomposed into plane-wave modes, and the three components with wavevectors 
that form an equilateral triangle may have different relative signs. The sign of the bispectrum tells you which combination 
of signs is more likely (on average gives a positive or negative product of the three modes) . A positive reduced bispectrum 
corresponds to being likely to have waves combining to have strong overdensities surrounded by larger areas of milder under- 
density. A negative equilateral bispectrum corresponds to being likely to have concentrated underdensities surrounded by areas 
of milder overdensity. Note that in 3D the figures extend into the page, and hence the positive bispectrum corresponds to 
concentrated overdense filaments surrounded by larger areas of milder underdensity. 




FIG. 2: A snapshot of non-linear large-scale structure from the millennium simulations [T5]. Dynamical non-linear collapse 
of very dense filaments (surrounded by milder underdensities, voids) generates a large positive roughly equilateral density 
bispectrum. 



during the growth of large-scale structure, as shown in the famous simulation of Fig. [2] Since it is the overdensities that 
are concentrated, not the underdensities, the non-linear large-scale structure density field will have a large positive 
equilateral component to its bispectrum (for a perturbation theory calculation see Ref. [IS])- 

Of course exactly equilateral triangles are a very special case, but there are many shapes that are close to equilateral 
and these will also look similar, but correspond to slightly elliptical concentrated overdensities or underdensities. 
As the bispectrum triangle becomes more flattened, these turn into a line, or in 3D concentrated overdensity or 
underdensity pancakes (planes); see Fig. [SJ Note that shapes that are qualitatively distinct in 3D may not be after 
projection into 2D: for example if an purely equilateral shape is present in 3D, projecting down to 2D will give 
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FIG. 3: As an approximately equilateral bispectrum triangle flattens, the round areas of overdensity become flattened into 
pancakes. In 3D a positive flattened bispectrum with fei = fc2 = ^3/2 corresponds to being likely to have overdense pancakes 
with larger mildly underdense planes in between. 

flattened contributions when the line of sight lies close to the plane of the triangle (slicing a 3D filament along its 
length gives a line of overdensity) . 

Since equilateral-form non-Gaussianity involves wavevectors of roughly the same magnitude, these modes would 
have left the horizon during inflation at roughly the same time. Non-linear dynamics prior to horizon exit during 
inflation can therefore generate a bispectrum with a significant equilateral bispectrum component. By an analogy 
with large-scale structure growth, one might imagine than any non-linear growth prior to horizon exit would require 
a low sound speed for the perturbations^; however in standard single-field inflation the rest-frame sound speed of 
the perturbations is exactly the speed of light. This prevents significant non-Gaussianity developing. However in 
extended models the speed of sound can be much lower, and in such cases significant equilateral non-Gaussianity can 
develop [Bj. Non-linear dynamics after horizon re-entry will of course also generate equilateral non-Gaussianity, for 
example second-order effects prior to recombination. Though small compared to the strongly non-linear growth of 
structure at lower redshifts, these signals might present an important source of possible confusion for small infiationary 
signals [14J. 

B. Squeezed triangles 

Squeezed triangles correspond to having one wavevector much shorter than the other two: in other words one 
large-scale mode and two much shorter-scale modes. The bispectrum is invariant under permutations of fci, A:2, ^3, so 
for squeezed triangles it is convenient to adopt the convention that we permute indices so that fci < fc2 < fca, and fci 
therefore always labels the large-scale mode. Sometimes people refer to the "squeezed-limit", meaning the limit as 
ki — > 0, but this is not really observationally relevant as very super-horizon modes are unobservable. By squeezed I 
will mean triangles with fc2, A:3 ^ fci, but the wavelength of fci not much larger than the horizon size today, so that 
the mode is still observationally relevant. 

First it is helpful to consider what a combination of two small-scale modes with k2 ^ — ks looks like: as shown in 
Fig. [4] the waves destructively interfere in some regions leaving little small-scale structure, but in other regions they 



^ The sound speed, given by c^ = 6p/Sp, measures the pressure perturbation Sp induced by a given density perturbation 5p. In structure 
growth pressures prevent gravitational collapse. The sound horizon roughly determines the Jeans scale below which pressure prevents 
collapse, so a low sound speed is required in order to have much growth on small sub-horizon scales. 
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FIG. 4; Two short-scale modes combine giving interference patterns. These look like a large-scale modulation in the amplitude 
of the small-scale modes, with the modulation having wavevector ki = — k2 — ka. 

reinforce each other giving a large small-scale signal. So this looks like a large-scale modulation in the small-scale 
power, where the wavevector of the modulation is given by — (k2 + ka). In a Gaussian field the signs of all the modes 
are independent, so of course there is no modulation on average. However if there is a non-zero squeezed bispectrum, 
there is a correlation between this modulation and the large-scale modes; for example see Fig. [51 A positive squeezed 
bispectrum means that where there's a large-scale overdensity there's likely to be more small-scale structure, and 
where there is an large-scale underdensity there is likely to be less small-scale structure. 



1. Angular dependence and the squeezed approximation 

So far we have described the bispectrum triangles in terms of the lengths of the sides. However for squeezed triangles 
in particular it can be useful to describe the triangles in a different way. Since the two small-scale modes are of similar 
wavelength, is natural to use a single number fc = |k2 — k3|/2 to quantify the scale of the small-scale modes [T^I16) . 
as shown in Fig. |6] The remaining free parameter can then be taken to be the angle (j) between the large-scale and 
small-scale wavevectors, so that a triangle is fully described by the three numbers ki,k,(f). Squeezed triangles have 
ki ^ k, and in particular I shall refer to the squeezed limit as having leading corrections of ©((fci/fc)^). Note that 
exchanging fc2 4-> k^ is equivalent to -(-)• tt — with the same k. 

It can also be useful to decompose the bispectrum depending on its angular dependence, i.e. writing^ 



6(fci,A:,<^)=^5^^, 



imcp 



(4) 



If the squeezed bispectrum is independent of (j), i.e. the orientation between the large and small-scale modes, it is 
called isotropic. In this case the bispectrum is fully determined by the m — component b^ j,. For 3D or statistically 
parity-invariant 2D fields than can be no sm{m(J3) dependence, and from rotational invariance the odd m components 
should vanish, but in general there can be an angular dependence cos(m(/)) with even m. This generally enters from 
ki •k2 and ki -ka dependence of the bispectrum, for example from gradient contractions, and hence angular dependence 



^ For bispectra in 3D, one could also expand into more directly orthogonal spherical harmonic Y^q modes; the argument is much the same 
so for simplicity I will stick with the 2D modes. 
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FIG. 5: Squeezed bispectrum: two small-scale modes with nearly-equal wavelength (^2 ~ fcs) interfere with each other, giving 
some regions with lots of small-scale power and others with destructive interference giving little small-scale power. The sign of 
the bispectrum tells you whether a region of high small-scale power is more likely to be associated with a large-scale overdensity 
or a large-scale underdensity. If the correlation is independent of the relative orientation (upper and lower figures have the 
same signal) the bispectrum is isotropic, but in general it is not. 




FIG. 6: A bispectrum triangle can be described by the lengths of the three sides, or alternatively by the length of the shortest 
side (corresponding to the large-scale mode), the length-scale of the long sides (short-scale modes) defined by fc = |k2 — k3|/2, 
and the angle (j> that measures the relative orientation of the long and short-scale modes. 



enters in the combination ki cos ((/>). A bispectrum expansion in ki/k therefore has angular dependence entering via 
powers of the dimensionless parameter e^ = (fci/fc) cos(0). Since any bispectrum must be symmetric under fc2 O k^, 
(and therefore •(-> vr — </)), only even powers of e^ can enter, and higher even powers of cos(0) (and hence higher 
even m) are suppressed by proportionally more powers of ki/k. Thus the squeezed expansion of a smooth scalar 
bispectrum is typically of the form 

1.2 

b{ki,k2, fcs) = A{ki, k) + [B{ki,k) + C(fci, k) cos(2(^)] -| + D(fci, fc, (l))0{k\/k^), (5) 

where A, B, C, D encode the scale-dependence of the particular physics involved. Hence unless there is very strong scale 
dependence the leading term is isotropic, and the next leading term only has isotropic and quadrupolar contributions. 
For some bispectra A = and the leading term then in general has both isotropic and quadrupolar components. 
Physically the reason higher angular dependence does not appear is because a small patch on a large-scale scalar 
modulating field will be accurately described by a field value and gradient, with the field value giving an isotropic 
change to the small-scale power, and the gradient defining the local basis for the quadrupolar dependence of the 
small-scale power. The leading term can only have anisotropic contributions if the modulation is not scalar. 

For example gravitational lensing of the CMB generates both m = part of the bispectrum (corresponding to 
large-scale lenses isotropically magnifying and de-magnifying the CMB), and also an to = 2 component from lensing 
shear (see e.g. Refs. [T51 [TTHTO] ). Primordial bispectra involving gravitational waves can also generate small m = 2 
components, since the modulation of the small-scale modes will depend on their orientation with respect to the 
anisotropic distortion produced by a large-scale tensor mode |2Qj . Similarly inflation models with vector fields can 
generate anisotropic bispectra (e.g. [21]). Purely scalar local modulations are expected to give isotropic bispectra, 
as we discuss further below, though leading terms from small gradients will give quadrupolar contributions (see e.g. 
Ref. Ea). 

The decomposition of the squeezed bispectrum into angular moments is conceptually useful for distinguishing 
different physical effects. In particular, in the squeezed limit bispectra with different to are orthogonal: an estimator 
for &™ J. should not be liable to confusion with a bispectrum 6™ j, if to, ^ m'. For example this is partly what 
allows CMB lensing to be easily distinguished from an isotropic local primordial bispectrum: the m = part due 
to magnification is a source of confusion, but the m — 2 lensing signal is distinctive and allows the lensing to 
be isolated and subtracted. The angular decomposition is also useful when considering secondary processing of 
primordial bispectra: if a statistically isotropic small-scale process affects the primordial modes it will not mix 
bispectra of different to. For example under gravitational lensing any primordial squeezed 6™ ;, is blurred out in k due 
to random small-scale lensing deflections, but the angular to dependence does not change [521 [H]. Since all bispectra 
are expected only to have the lowest m components in the squeezed limit, this also suggests that a general modal 
bispectrum decomposition can efficiently capture the angular dependence with only two angular modes for squeezed 
shapes, and since corrections are quadratic triangles do not need to be very squeezed for the squeezed approximation 
to be quite accurate. 

Models in which statistical isotropy is broken can also generate anisotropic bispectra in the more general sense that 
the bispectrum signal can then depend on the orientation of the triangle [551 [2S]j I '^iU not consider this possibility 
further here. 

2. General form of the squeezed bispectrum 

Consider the case where the field being observed T can be calculated from some non-linear function of a set of purely 
Gaussian fields {Xi} that we can write as a vector X. For example the Gaussian fields might be combinations of 
small scale and large-scale perturbation modes, different inflation field perturbations, lensing potentials, linear matter 
densities, velocity fields giving redshift distortions, etc. Since we are assuming underlying Gaussian fields, they are 
fully described by their covariance C (a matrix of power spectra) . If we are interested in the bispectrum of a field 
T{X) (where the tilde denotes that it is non- linear), it can therefore be calculated by integrating out the Gaussian 
fields 

(f(ki)T(k2)f(k3)) = |dX^^^^H^|^^f(ki)r(k2)f(k3). (6) 



^ In 2D, see Ref. |23| . For bispectra in 3D, one can expand in spherical harmonics Yi^q, and the modes with different L will then be 
orthogonal. 



For a squeezed bispectrum the large-scale field r(ki) is at a much larger scale than the two small-scale fields T(k2), 
T(k3), and hence is often well approximated as being Gaussian (a linear function of the underlying Gaussian fields). 
For example in large-scale structure the large-scale mode in the bispectrum will be nearly linear if ki < 0.05Mpc~^, 
similarly in the CMB large scales are accurately linear but smaller scales have significant non-linearities due to lensing 
and other effects. The approximation T(ki) « T'(ki) is the linear short-leg approximation [15 , which is very accurate 
in some cases. Writing T(ki) as a linear combination of Xi(ki) in Eq. (|6|, T(ki) = MiXi(ki), and then writing 
X,;(ki) as a functional derivative of the exponent, 



(r(ki)r(k2)f(k3)) ^ |dx '^P ^|2!c|V2 '^^ M.^.(ki)f (k,)f(k 



3) 



Then since Cij — PxiX , we have MiCij — Ptx and hence 

(f(ki)f(k2)f(k3)) « (r(ki)f(k2)f(k3)) - PTxSkl) {jxj^ (r(k2)f(k3))y (8) 

Thus the bispectrum depends on the correlation of the large-scale field with the modulating fields PTXi{ki), and is 
proportional to the response of the small-scale non-linear modes to changes in the large-scale modulation. The result 
of Eq. ([8]) is fully non-perturbative and only relies on the linear short-leg approximation, not an extreme squeezed limit 
(if necessary it can be generalized out of the linear short-leg approximation by including higher derivative terms) . In 
the case of CMB lensing it is possible to calculate the response term essentially exactly non-perturbatively [TS1 ; more 
generally as ki/k — >■ 0, fc2 — > ^3, and the response term just describes how the small-scale power spectrum changes 
with a different large-scale background modulation. 

3. Local scalar modulations 

Perhaps the simplest way to generate a bispectrum is from a purely local modulation of a Gaussian field. For 
example 

x(x)=xo(x)[l + V(x)] (9) 

where both xo and ^ are Gaussian fields. I will assume that \ip\ <C 1 so that the modulation is small, consistent 
in the primordial context with the relatively tight observational bounds on non-Gaussianity of the CMB. To have 
a non-zero bispectrum the fields must be correlated, so Pxo4> t^ 0- The correlated part of ^ is then conventionally 
written as I/nlXOi where /nl is an amplitude determined by the degree of correlation (and hence governing the size 
of the bispectrum). It is then straightforward to calculate the bispectrum to leading order in the modulation, giving 

b{k,,k2,ks) = 2^/NL[P(fcl)P(fc2) + P{kl)P{k3) + P{k2)P{k3)] (10) 

5 

where P{k) is the power spectrum of xo- 

For primordial perturbations usually nearly scale- invariant fields are of interest, what does that mean in terms of 
/c-dependence of P{k)l The total variance at any point in an N-dimensional field is ex J d^kP{k) ex J k'^ P{k)dlnk. 
A scale-invariant spectrum has equal contributions to the variance from each dimensionless d In A: interval and hence 
has P{k) ex k~^ . If the fields are close to scale invariant, for squeezed triangles with ^2,^3 ^ ^i it follows that 
P{k2)P{k^) ^ P{ki)P{k2) ~ P{ki)P{k-i). In fact doing an expansion in ki/k for a scale-invariant spectrum we have 

b{ki,k2,k3) « 2^/NLP(fcl)(P(fc2)+P(fc3)) 

5 



= ^lhLP{ki)P{k) 



kiV 9 + 15cos(20) 
TJ 16 "*" 



(11) 



This demonstrates that, in the squeezed limit, the local modulation gives an isotropic bispectrum as expected. The 
leading angular dependence is quadrupolar (with m = 2) and only enters at 0{{ki/k)'^), corresponding to anisotropy 



9 

generated by large-scale gradients in the modulating field: to leading order all the small-scale modes have their 
amplitude modulated the same way, independent for their orientation. This remains true for smooth nearly scale- 



invariant spectra, in general Eq. (11 1 just picks up contributions to the 0{{ki/k)^) term that depend on the derivative 
of the power spectrum. 

Local modulations appear in large-scale structure, where the modulating field is a large-scale density mode. Above 
linear order this will affect how small-scale perturbations grow, and hence modulate their amplitude. However the 
matter power spectrum on scales larger than the horizon size at matter radiation equality (the 'turnover') has rapidly 
decreasing power, so for the density field this effect is only important for squeezed triangles where all scales are smaller 
than the turnover (in the fci — )■ limit the density bispectrum is zero, see e.g. Ref. [in])- On the other hand if there 
is primordial non-Gaussianity in the form of a local modulation by a gravitational potential, this will also modulate 
the amplitude of the small-scale power, and hence modulate the number of objects that reach the critical density 
to form a galaxy. Since the potential is scale invariant on large-scales (proportional to the comoving density over 
k'^ by the Poisson equation), the large-scale potential modulation of the small-scale density leads to an observable 
large-scale modulation in the number density of galaxies that form. Exactly how this shows up in the power spectrum 
and higher-point functions is a little subtle, see Refs. [TUJ [571 HH] for details. The simplest effect is actually in the 
galaxy power spectrum: local primordial non-Gaussianity leads to a scale-dependent bias on large-scales. 

Local modulations are also potentially important in the early universe: once large-scale fluctuations are well outside 
the horizon they can give a somewhat different local evolution in different Hubble patches, which can locally modulate 
what happens to the smaller-scale modes. 

4- Squeezed bispectrum from inflation 

For squeezed bispectra, the small-scale modes leave the horizon during inflation significantly later than the large- 
scale modes. The large-scale modes therefore effectively modify the background seen by the small-scale modes as 
they leave the horizon and beyond. In single-field inflation there is a one-to-one mapping between the background 
Hubble parameter and field values, so the conditions on the surface where a small-scale mode leaves the horizon at 
^phys — H are locally identical to what they would be without the large-scale mode [ini HO] • However when we observe 
angular scales on the sky, we can look at many different Hubble patches; since these patches all have slightly different 
large-scale field values they have undergone slightly different expansion histories: there is a local perturbation to the 
scale factor in each patch. 

Is this observable even though the physics at horizon exit at all points is the same? It might be helpful to consider 
the analogy of the CMB fiuctuations: if we approximate recombination as sharp so that the universe goes immediately 
from opaque to transparent, the CMB photons we observe are coming from an equal temperature surface (given by 
the critical temperature at which electrons and protons recombine). However we do not see the same temperature! 
On small scales this is because in locations where there is an overdensity — which is hotter — the universe locally 
has to expand for longer to reach the recombination temperature. The photon we observe has therefore seen less 
expansion (and hence redshifting) between recombination and us, and hence appears hotter. In this picture it is the 
local scale-factor perturbations at recombination that give rise to the small-scale observed temperature anisotropy, 
even though the surface itself is equal temperature. In the case of the inflation fluctuations the picture is similar: the 
large-scale modes affect the amount of expansion between horizon exit and observation, and hence the same physical 
scale has been stretched by slightly different amounts in different places. The difference in size between the many 
different Hubble patches is in principle observable; an observed angular scale in different directions corresponds to a 
different physical size exiting the horizon, and hence a slightly different power spectrum amplitude. 

However the effect is small, most obviously because the fluctuation amplitude is 0(10"^), so the modulation is not 
expected to be at larger than 0(10"^), at least if the inflaton potential is smooth. The change in power with scale 
is also small because the spectrum is expected to be nearly scale-invariant: from Eq. ([8]), for a power law spectrum 
k^P{k) ex fc"""^ we have the well-known result for the curvature perturbation bispectrum in the squeezed limit |20j 
(see the Appendix for derivation) 

5(fci, fc2, fca) = fe^ifc « -^mo(n. - l)F(fci)F(fc). (12) 

For rig ^ 1 this corresponds to |/nl| <C 1, so a detection of primordial squeezed non-Gaussianity at a level higher 
that this would rule out essentially all single-fleld inflation models. However in models with multiple fields the super- 
horizon evolution can be modulated in more significant ways by auxiliary fields, and larger non-Gaussian signals are 
possible (for a review see Ref. [5]). 
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FIG. 7: A trispectrum shape is defined by four wavenumbers ki..k4, and the 'diagonals' Ki2,K2s- Lines with these lengths 
form the edges of a tetrahedron. 

IV. TRISPECTRUM 

The trispectrum is the harmonic four-point function. As with the bispectrum, statistical isotropy and homogeneity 
means that the four wavevectors must sum to zero, and for a parity invariant ensemble the shape is fully determined 
by the lengths of six lines joining the four points. Gaussian fields have a trivial four-point function given by products 
of the power spectrum, so the trispectrum T is defined as the more interesting connected part: 

(T(ki)T(k2)r(k3)r(k4))c « S{k, + k2 + ka + k4)r(fci, fc2, h, k^- K.^.K^s). (13) 

Here ki..k4 are the lengths of the wavevectors, and K12, K23 are the lengths of the two 'diagonals'; see Fig. [7| Note 
than in 3D the points do not all have to lie on a plane: in general the 6 connecting lines define the edges of tetrahedron. 

A. Tetrahedral trispectrum 

The most symmetric case is when the tetrahedron defined by the connecting wavenumbers is regular; a field with 
this kind of 'tetrahedral' trispectrum non-Gaussianity is shown in Fig. [8] The positive sign corresponds to having 
concentrated regions of overdensity and underdensity, with weaker perturbations in between. The histograms show 
the point kurtosis, with sign depending on the sign of the trispectrum. Such a trispectrum could be generated by 
local dynamics, with strongly non-linear local processes potentially giving a large signal, or inflation models with low 
sound speed. But note that because the trispectrum is invariant under T -^ — T it does not describe skew, and hence 
on its own does not describe the effects of gravitational collapse (which generates a large equilateral bispectrum) . Of 
course the tetrahedron does not need to be regular, and non-linear dynamics will generally also generate a tetrahedral 
signal with similar but not equal sides, corresponding to more stretched-out regions of overdensity and underdensity. 

In 2D (e.g. on CMB sky) the tetrahedron becomes a quadrilateral. Specific shapes of interest include kites and 
trapeziums^ |29j , symmetric examples of which are shown in Fig. pi Exact parallelograms are useful for calculating 
the non-Gaussian contribution to the covariance of the power spectrum, but near-parallelograms are equivalent to the 
squeezed trispectrum that I discuss in more detail below. 



* Note that I use 'trapezium' in the international (non-US) English sense; a quadrilateral with two parallel sides, with adjectival form 
'trapezoidal'. 
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FIG. 8: Density contours of a three-dimensional field with 'tetrahedral' trispectrum, corresponding to the sum of four modes 
of equal wavenumber, with wavevectors determined by the directions of the vertices of a regular tetrahedron. The figures on 
the left are for positive trispectrum (positive kurtosis as shown in the histogram of point values); on the right for negative 
trispectrum (with negative kurtosis). A positive tetrahedral density trispectrum corresponds to having localized regions of 
concentrated overdensity and concentrated underdensity surrounded by larger volumes that are only weakly perturbed. 
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FIG. 9: Two dimensional fields with trapezoidal and kite-like trispectra of both possible signs. The specific shape trapezium 
shape has k2 = k^ — k-^ = fei/2, and the kite is a cyclic quadrilateral with 2fci — 2fc2 = fcs = fc4. 
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FIG. 10: Positive trispectra with one short diagonal (diagonal squeezed), corresponding to a large-scale modulation of small- 
scale power. In 3D fields the wave vectors can also go into the page 



B. Squeezed-diagonal trispectrum 



As in the case of the bispectrum, the squeezed limit is of particular interest. The diagonal-squeezed trispectrum 
corresponds to one of the diagonals of the tetrahedron being much smaller than the other sides, e.g. i^ ^ fci, ^2, fcs, fc4. 
The tetrahedron can be considered as two triangles having side K being stuck together; hence large-scale wave number 
K is the wavenumber of the beats seen in Figj4j and corresponds to the wavelength of a large-scale modulation in the 
small scale power. This is illustrated in Fig. 10 for wavevectors lying in the same plane; the shape of the tetrahedron 
determines the relative alignment of the modulation and the small-scale modes. In general the squeezed shape depends 
on both the diagonals, but in the two-dimensional case, where all the wavevectors lie in a plane, the squeezed-diagonal 
trispectrum only depends on the short diagonal K since the other diagonal is then determined from the lengths of 
the sides. Note that an arrangement of wavevectors into an approximately planar parallelogram is equivalent to this 
configuration, since the wavevectors can be re-ordered so that the nearly parallel vectors are consecutive. 

I'll discuss the general form of the theoretical prediction for a squeezed-diagonal trispectrum later in Sec. [Vj The 
local modulation of Eq. Q can produce it, with all small-scale modes being modulated in the same way. Such a 
trispectrum is positive (because the modulation field always has positive power spectrum) and isotropic: there is no 
dependence on the relative alignment of the modes, and the signal in all the configurations shown in Fig. [lO] should 
be the same. In general however there can be angular dependence; for example in the 2D case of the CMB (where 
all wavevectors must lie in a plane), CMB lensing shear generates an anisotropic squeezed trispectrum [3D]. Similarly 
large-scale gravitational waves produce small orientation-dependent distortions in the small-scale modes [SU [32] . 

Note that the modulation field here is only observed via its effect on the small-scale modes; as such only the 
amplitude of the modulation matters, not its correlation to the large-scale field. Further note that if there is a local 
modulation giving the bispectrum, it must also give a squeezed trispectrum: if there's a spatial modulation of small- 
scale power correlated to the large-scale temperature giving a bispectrum, the power of the in large-scale modulation 
field gives a trispectrum. The correlation of the modulation can be between —1 and 1, which limits the size of the 
bispectrum compared to the trispectrum, with the bispectrum signal being maximized when the modulation field 
is totally correlated to the large-scale field. Conventionally the size of any primordial isotropic squeezed-diagonal 
trispectrum is parameterized by a parameter tnl and there is then an inequality tnl ^ (6/nl/5)^ in the squeezed 
limit 33 . If the amplitude of the modulation depends on K, there is also an equivalent relation between the bispectrum 
with ki = K and the trispectrum at each K, which just says that the correlation of a modulation at any scale with 
the equivalent large-scale field mode must be less than one. 

The diagonal-squeezed trispectrum is hard to measure. Firstly, because if it is due to an 0{e) modulation, it is 
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One-leg squeezed trispectra |fei| « |/(:2l~l'i^3h|fc4l 

Trispectrum>0 Trispectrum<0 

• •••••••••••♦••••••• 

• .•....••••••••...... M.M.l.M.M.X.I.l.l.l.l.l.M.l.M 



FIG. 11: Trispectrum with one short wavenumber, fci ^ ^2,^3,^4, corresponding to a large-scale modulation of a small-scale 
nearly-equilateral bispectrum. In 3D fields the modulation can also go out of the plane so the sign of the density changes along 
each filament. 

proportional to the power in the modulation which is a tiny 0{e^), so a larger modulation is required to have an 
observable signal. Secondly because of cosmic variance: the finite number of small-scale modes available with finite 
resolution observations. If we imagine that the field consists only of large-scale and small-scale modes, to measure 
the modulation of the small-scale power we need to divide up the space and measure the small-scale power as a 
function of position. Each box in space will have a sample variance error oc 1/n in the small-scale power, where 
n is the number density of the small-scale modes. For boxes large compared to the small-scale modes this error is 
independent between boxes, and the sample variance therefore looks like white noise. This translates into a white 
cosmic variance noise on the estimate of the modulation field. If the modulation field is nearly scale-invariant, 
the sample variance noise therefore dominates on small-scales. In that case only the very large-scale modes of the 
modulation field have any chance to be measured, and in fact the signal to noise (for small signals) is dominated by 
the very largest modulations (in the CMB, almost all the signal is in modulations with L < 10 J4J). By contrast 
the squeezed bispectrum is the correlation of the modulation with the large-scale field, where the large-scale field is 
typically measured very accurately; so strongly correlated modulations are generally much more easily seen in the 
bispectrum than trispectrum. However in cases like CMB lensing, where the correlation drops rapidly with scale and 
the convergence power spectrum is very blue, the trispectrum is easier to detect. Note that if the modulation field 
starts to be well measured, increasing the number of small-scale modes to lower the sample variance noise does not 
help any more: the accuracy of the bispectrum and trispectrum measurement is then limited by the cosmic variance 
of the large-scale modulation modes [U [TS] . 



C. One-leg and two-leg squeezed trispectra 



For the trispectrum there is another squeezed limit to consider: when one of the edges in the tetrahedron is 
much shorter than the others. A field with this kind of trispectrum is shown in Fig. |11[ The large-scale mode with 
ki <^ ^2,^3, ki modulates the small-scale modes, and if ^2 ^ fcs ^ k^ this look like a modulation of a nearly equilateral 
bispectrum. As shown in Fig. 11 the small-scales look like concentrated overdensities (positive equilateral bispectrum) 



or concentrated underdensities (negative equilateral bispectrum), where the sign of the small-scale bispectrum is 
correlated to the large-scale modulating field value. A positive one-leg squeezed trispectrum corresponds to regions 
with small-scale concentrated overdensities being more likely where there is a large-scale overdensity, and concentrated 
underdensities being more likely where there is a large-scale underdensity. For negative trispectrum the association 
is the other way round. 

If the bispectrum that is being modulated is local rather than equilateral, the signal will peak when ^2 ^ fcs, k^ and 
hence in a two- leg squeezed trispectrum shape with ki,k2 <S^ k^jk^. A trispectrum of positive sign can be obtained 
from a correlated linear local modulation as in Eq. (|9]), and hence is produced by local primordial non-Gaussianity (in 
addition to a diagonal-squeezed trispectrum). However the negative sign cannot be produced in this way. However a 
local quadratic modulation can give either sign, for example 



x = xo(i + gNLX§) 



(14) 
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FIG. 12: Examples of two dimensional fields with flattened trapezoidal trispectrum. In the top figure the waves have equal 
amplitude, and in the flat limit where ^2 = ^3 = ^4 = ki/3 the negative sign trispectrum would come from the first two terms 
in the harmonic expansion of a square wave. The lower figure has wavenumbers in the ratio 1, 3, 5, 7 with equal amplitude; if 
the amplitudes were in the ratio 1, 1/3, 1/5, 1/7 the positive trispectrum would come from the first four terms in the harmonic 
square-wave expansion. Cosmic strings are expected to produce the top trispectrum with negative sign, and the bottom one 
with positive sign (in agreement with Ref. |29]), corresponding to having sharp real-space line-like changes in temperature. 

where xo is Gaussian and (7nl quantifies the size and sign of the modulation. This form of modulation does not produce 
a bispectrum, but can produce the two-leg squeezed trispectrum of either sign, which can be verified by explicitly 
calculating terms in the four-point correlation. A negative two-leg squeezed signal is therefore clear signature of 
this kind of primordial non-Gaussianity. Unfortunately it is rather hard to measure: as we've seen, the signal is a 
correlation of a bispectrum modulation, so estimating the amplitude effectively requires estimating the bispectrum in 
different local patches and then correlating its spatial dependence with the large-scale field. The cosmic- variance noise 
is therefore much larger than when estimating a constant primordial equilateral bispectrum and current constraints 
are rather weak |35| . 

The general form of the theory one-leg squeezed and two-leg squeezed trispectra follows an analogous argument to 
that leading to Eq. (Is]): in the linear one short-leg approximation: 



(f(ki)f(k2)f(k3)f(k4)) « (r(ki)f(k2)f(k3)f(k3)) =PTX,(fcl) 



6x,{k,y 



(f(k2)f(k3)f(k4) 



(15) 



The single-field inflationary prediction is very small because this trispectrum depends on changes in the bispectrum, 
and the bispectrum itself is very small |36j . The signal in large-scale structure is expected to be much larger and of 
positive sign, since non-linear growth gives an equilateral bispectrum, and the amount of non-linear growth will be 
modulated by large-scale densities. The two-leg squeezed trispectrum can be calculated similarly assuming linearity 
in the two short-legs, and can be useful for example when calculating analytic approximations for the ^nl trispectrum 
expected in the CMB ^5] . 



D. Flattened quadrilateral trispectrum 



We have already discussed trispectra with a squeezed diagonal, which if all the wavevectors are nearly aligned is an 
example of a flattened quadrilateral (e.g. the lower flelds in Fig. 10). In 2D the flattened limit corresponds to having 



a real-space symmetry along lines. Other cases of interest are the flattened trapezoidal shapes shown in Fig. 12 
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If there are cosmic strings, the Kaiser-Stebbins effect [37] can produce hne-hke discontinuities in the CMB temper- 
ature. This is because the strings are very thin, but moving fast: in their rest frame they see a large CMB dipole, 
which is lensed by the string stress-energy, giving a change in observed temperature across the string. So strings 
produce sharp Hues in the CMB where the temperature suddenly changes from hotter to colder. There is therefore 
power on very small scales (big change in temperature across the string), with little power in between strings; this is 
like a very localized modulation of small-scale power, corresponding to a squeezed-diagonal trispectrum shape. 

An array of strings with appropriate velocities would give a signal in the CMB that looks like a plane square wave, 
with the line-like discontinuities corresponding to the edges of the steps. Expanded into Fourier modes, a square wave 
with wavevector in the x direction gives 

y cos^j2ri + l]kxl ^ ^Qg(^^) _ 1 cos(3A:a;) + ^ cos(5fca;) - ^ cos(7fca;) . . . , (16) 

■^-^ 2n + 1 3 5 7 

n 

SO there is a relation between the amplitudes seen in the modes with wavenumber fc, 3fc, 5k and 7k. For example 
T{k) = -3T(3fc), and thus we expect {T {k)T {3k)T {3k)T {3k)) < 0; likewise we expect {T{k)T{3k)T{5k)T{7k)) > 0. 
This can be seen in Fig. |12[ where only one of the signs in each configuration looks roughly like an array of line 
temperature discontinuities as expected from cosmic strings. So cosmic string networks will produce negative trispectra 
of some flattened shapes, along with trispectra in other flattened (and less-flattened) configurations with positive 
sign [29l [33 • Note that cosmic strings produce a small bispectrum (but significant trispectrum) since the velocities of 
long strings producing the Kaiser-Stebbins signal can be in any direction, so that to a first approximation T(x) and 
— T(x) are equally likely. 

Other flattened quadrilateral shapes will of course be produced by other ID density functions. For example delta- 
function line sources contribute equally to all wavenumbers in the preferred direction, so there would be a positive 
flattened trispectrum with wavenumbers in the ratio 1:2:3:4 where k4 + ki = — k2 — ks. 

V. STATISTICAL ANISOTROPY AND MODULATION RECONSTRUCTION 

We've seen that large-scale modulating fields can give rise to a squeezed bispectrum and trispectrum when we 
average over all possible field values. However we could ask a slightly different question: in our universe there must be 
a particular realization of the large-scale modulation field, can we learn what it is? Looking at Figs. [5] and [TO] how to 
do this seems obvious: the small-scale power is modulated, so if we estimate the power spectrum in different places, 
the spatial dependence of the power spectrum should trace out the modulation. 

Imagine that there are modulation scalar modes X that we want to learn about. If they are large-scale modes, the 
analog of the linear short-leg approximation here is that products of small-scale fields only depend linearly on X (but 
possibly non- linearly on all the other fields^). For a pair of small-scale modes we therefore have 

f(k2)f(k3)«f(k2)f(k3)U=0 + y'dKX(K)*^^^(f(k2)f(k3)). (17) 

If we then average over the conditional distribution for the fields given fixed X, for k2 7^ — ka we have 

(f (k2)f (k3))p(t|^) « J dKX(K)* ^^^^ (f(k2)f (kg)) ^ , (18) 

where for a result to linear order in X(K) the expectation on the right hand side can be evaluated averaging over 
all the fields. If we correlate this with T(ki) and average over X this recovers the general form of the squeezed 
bispectrum given in Eq. (|8|. From statistical homogeneity the wavevectors in the expectation value over all fields 
must form a triangle, so this term is proportional to a delta-function, and we can define 



6X{K)* 



r(k2)T(k3) ) = A{K, fc2, k3)SiK + k2 + fc3), (19) 



where A{K, ^2,^3) encodes how the small-scale modes change with the large-scale modulation. Hence averaging over 
all the other modes recovers something just proportional to the modulation ^(K) [ISt i39i : 

(f(k2)f(k3))p(f|;,) «^(i^,fc2,fc3) ^(K)*|K=-k.-k3 ■ (20) 



Note that the wavenumber dependence of X can be restricted as much as we like, so this only requires linearity in a single mode: on 
the CMB the dependence can be non-linear in all modes as long as it is linear in the single large-scale l,m mode being considered. 
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This tells us that by averaging over a quadratic combination of all the observed small-scale modes (which all see the 
same large-scale modulation realization), we can construct an estimator X of the modulation field X. This is the more 
formal statement of the obvious idea that measuring the small-scale power as a function of position should trace out 
the large-scale modulation. The use of such quadratic estimators has proved useful in CMB studies for reconstructing 
fields in wide classes of statistically anisotropic models [301 EI] , and for lensing reconstruction (where the modulation 
field is the large-scale lensing potential |42l |43]; for reviews see Refs. [44l |45]). Of course in observations there are 
only a finite number of small-scale modes to average over, and this leads to a cosmic variance reconstruction 'noise': 
random fluctuations in an isotropic field will look anisotropic, which is hard to distinguish from a modulation-induced 
anisotropy with only a small number of observed modes. The CMB lensing signal is large enough to be relatively easily 
detected (as seen recently by ACT [46]; see [43] for other forecasts), but primordial signals are much more difficult 



since only the very largest-scale modulations are not swamped by reconstruction noise (as discussed in Sec. IV B). 

For a fixed modulation the field can look anisotropic and inhomogeneous. However if the modulation is a statistically 
homogeneous and isotropic field, after averaging over the modulation everything is again statistically homogeneous and 
isotropic. The bispectrum and trispectrum then respectively quantify the correlation and variance the modulation. 
Specifically, correlating the modulation reconstruction with r(ki) yields an estimate of the bispectrum [HHTTIHS] : 

(r(ki)f(k2)f(k3)) = (r(ki)(f(k2)f(k3))p(^i;,)^^^ 

= S{ki+k2+k3)PTx{ki)A{h,k2,h), (21) 

which will be non-zero if the modulation and large-scale fields are correlated, Ptx y^ 0. If the field for fixed modulation 
is Gaussian, the power spectrum of the modulation is proportional to the squeezed-diagonal trispectrum since 



(r(ki)T(k2)T(k3)T(k4)) = ^(T(ki)T(k2))p(^|;,)(r(k3)r(k4))p(^l^) +perms)^ 

« ,5(ki + k2 + k3 + k^)Pxx{K)A{K, ki,k2)A{K, h, k^) (22) 

where in the last line I assumed K — |ki + k2| = |k3 -I- k4| ^ ki,k2,k^,ki is the short side of the tetrahedron (but 
K ^Q) and a nearly scale-invariant modulation so that Pxx{K) ^ Pxx(|ki + k3|),Pxx(|ki + k4|). 

For squeezed non-Gaussianity where Eqs. (21 1,(22) hold, since P|-y < PxxPtt there is an inequality between 



the amplitude of the trispectrum with ki,k2,k^,k4,K and the amplitude of bispectra with wavenumbers K,ki,k2 
and Kjks^ki (and the power spectrum Ptt)'- if you have a squeezed bispectrum you ine vitabl y have a non-zero 



trispectrum. For the local model this gives the tnl > (6/nl/5)^ inequality discussed in Sec. IV B (for full derivation 
in the case of scalar modulations see Ref. [49], though an equivalent inequality holds more generally). In terms of 
optimal estimators for the modulation (summing over all the contributing wavenumbers), this inequality is often 
obvious; for example for CMB lensing it just states that (C; )^ < Cf^Ci , where ip is the reconstructed lensing 
potential and T is the CMB temperature. There are similar inequalities involving higher-point functions; for example, 
a one-leg squeezed trispectrum corresponds to a correlated bispectrum modulation, so the bispectrum modulation 
must have some variance and hence there must be a non-zero 6-point function. 

VI. RANDOMLY LOCATED FEATURES 

If there is some particular fixed feature that can randomly appear on the sky, if the location is unknown, then the 
location-marginalized model is still statistically isotropic — there is no preferred direction or orientation of the feature. 
As such, the feature (or features) will give a signal in non-Gaussianity searches. For a simple example, say there can 
be a circularly-symmetric feature on the sky with some radial profile T{r) (for example as might be produced by a 
texture [SD] or inflationary bubble collision [51] [55] , in addition to other signals) . If it is centred at location c then 

T(k) = / |^T(x - c)e-*'^-^ = e-*"^-" f !:^:^r(j.)e-»fe'-™se ^ ^-^k■c f rdrT{r)Jo{kr) 

= e-'^-%k) (23) 

where t{k) is the Bessel-transform of T{r) defined by the previous line. Then if we integrate out the random location 
of the centre the n-point function is 

(T(ki)...r(k„))oc /dcJ]t(A:,)exp(-ik, •c)(x5|^k,J J]i(fc,). (24) 
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This is a particularly simple separable form of n-point function. As a specific simple case, if the feature is a constant- 
temperature disk of radius R, then 

(T(k,)...T(k„))cx<5(^^k.^n^^- (25) 

If the template is not rotationally invariant, then of course the orientation also has to be integrated out. If the feature 
changes sign under a rotation this leads to all of the odd n-point functions being zero, as for example in the simplified 
example of line-like temperature changes due to cosmic strings discussed in Sec. |IVD[ If the amplitude of the feature 
is also random and symmetric about zero, integrating out the amplitude will also give zero for all the odd n-point 
functions. 

If there is a feature with random sign, but there is only one on our sky, on our sky it will have some particular 
sign, and hence contribute to the bispectrum estimator in that realization. Is this consistent with prediction for 
the bispectrum being zero? Actually yes, because the variance of the bispectrum estimator is also sensitive to non- 
Gaussianity, and in this case the apparent signal in the bispectrum estimator can be interpreted as the variance of 
the bispectrum being inconsistent with the Gaussian expectations (in other words the large bispectrum estimator in 
our one realization is a detection of the six-point function) . In such cases it may be more natural to think about the 
conditional distribution, where the sky is anisotropic for a given fixed feature location. 

Localized features are of course unlikely to be best detected by looking at only the bispectrum and trispectrum, 
but if they are present then their contributions to these n-point functions can be important, and hence must be 
distinguished from other signals of interest. 

VII. CONCLUSIONS 

Non-Gaussianity studies can extract significantly more information from cosmological data than just the power 
spectrum. I've shown what the various types of qualitatively different signal look like in real space, which can be an 
intuitive aid to understanding the origin of possible non-Gaussian signals and which sources of non-Gaussianity might 
be confused with each other. I gave some general results for the form of the squeezed bispectrum and trispectrum 
and the relation between them. Squeezed configurations are especially interesting configuration to study in terms of 
learning about the physics of inflation. The angular dependence of any the squeezed signal can be used to identify 
different physical effects. Local modulations as produced by multi-field inflation are isotropic, but anisotropic models, 
or late-time effects like CMB lensing, can generate a distinctive anisotropic squeezed bispectrum. It can also be 
useful to think about non-Gaussianities in terms of statistically anisotropy; for a fixed modulation field the universe 
looks anisotropic, and the large-scale modes of the modulating field can be reconstructed by using observations of 
many modulated small-scale modes. However reconstructed modulations, and non-Gaussianity estimates generally, 
are often limited by cosmic variance: the finite number of modes that are available to distinguish random fluctuations 
in a realization of a purely Gaussian held from genuine non-Gaussianities. This is usually the limiting factor when 
trying to detect small primordial signals. 
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Appendix A: Squeezed bispectrum from single-field inflation 

For completeness I show how the well-known result for the single-fleld squeezed-limit bispectrum [THl [Ml IS2 can 
be obtained from the general form the squeezed bispectrum given in Eq. (pi) . In the case of single-fleld inflation the 
modulating field is the large-scale curvature perturbation C(ki). If fci ^ ^2, ^3 the background scalar metric seen by 
the small-scale modes can be written 

ds^ = dt^ - e^^a{tfdTi^ = dt^ - a{tfdx^ (Al) 
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Here ( is essentially a local scale-factor perturbation, and in the second step we took the squeezed limit so that C(x) = ( 
is constant over the region of interest, and x = e''X is a re-scaled coordinate. In terms of x the metric just looks like 
a normal background FRW universe, so the non-linear small-scale mode can be obtained from C(x) — C(^) — C('3^x) 
where ({x) is the usual linear fluctuation obtained in single-field inflation. The response of the small-scale mode to 
the large-scale perturbation is therefore 
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It is possible to proceed non-perturbatively, but for small inflation perturbations it is not necessary, so we can evaluate 
the response to zeroth order in C giving 
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We then have 
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1 



fco d In h 



-iklPcdk2)) 



(A6) 
(A7) 

(A8) 



For a power law spectrum this gives Eq. (12 1 in the main text. Ref. 
gradients in the large-scale field. 
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